Charge transport through bio-molecular wires in a solvent: Bridging molecular 
dynamics and model Hamiltonian approaches 
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We present a hybrid method based on a combination of quantum/classical molecular dynamics 
(MD) simulations and a model Hamiltonian approach to describe charge transport through bio- 
molecular wires with variable lengths in presence of a solvent. The core of our approach consists 
in a mapping of the bio-molecular electronic structure, as obtained from density-functional based 
tight-binding calculations of molecular structures along MD trajectories, onto a low dimensional 
model Hamiltonian including the coupling to a dissipative bosonic environment. The latter encodes 
fluctuation effects arising from the solvent and from the molecular conformational dynamics. We 
apply this approach to the case of pG-pC and pA-pT DNA oligomers as paradigmatic cases and 
show that the DNA conformational fluctuations are essential in determining and supporting charge 
transport. 
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Can a DNA molecular wire support an electrical cur- 
rent? The variety of partially contradictory experimen- 
tal results obtained in the past years P, 0,M 0, 0, H, 0| 
has hinted not only at the difficulties encountered to 
carry out well-controlled transport measurements, but 
also at the strong sensitivity of charge migration to 
intrinsic (base-pair sequence, internal vibrations) and 
extrinsic (solvent fluctuations, molecule-electrode con- 
tact) factors. Recently 0, 0|, two experimental groups 
have measured similar high electrical currents on the 
order of 50-100 nA despite the fact that the electri- 
cally probed base sequences and lengths were rather 
different. In spite of considerable theoretical research, 
the dominant mechanisms for charge transport through 
DNA wires have not been, however, fully elucidated, 
see e.g., Ref. I sj fo r a recent review. Electron transfer 
experiments 9| , Il0| a s well as related theoretical stud- 



ies [nl [H, [Tj, 111, 111, [B, [l3, [M, 19, 20, 21, 22., m, ill 
have clearly pointed out the crucial role of dynamical 
fluctuations in favouring or hindering hole transfer. We 
may thus expect that this may also be the case for charge 
transport. Studies based on model Hamiltonian formula- 
tions describing disordered sequences 1251. |26| or the cou- 
pling to dynamical degrees of freedom |l8l.ll9l. 20l . 21, 2^ 
involve many free parameters which are in general diffi- 
cult to determine for realistic situations. First-principle 
calculations 27, m, m, H performed on static 

structures provide, on the other hand, orders of mag- 
nitude for the electronic coupling parameters but can 
hardly deal with the coupling to dynamical degrees of 
freedom. The inclusion of dynamical effects in quantum 
transport calculations has only been addressed in few es- 
pecial cases [l3, [16, 33, 34] in a systematic way. Thus, 
a general approach able to combine dynamical informa- 



tion drawn from a realistic description of bio-molecular 
conformational dynamics with a treatment of quantum 
transport is highly desirable. 

In this Letter, we present a study of charge transport 
through bio-molecular wires with different lengths by us- 
ing a hybrid approach based on a mapping of the time- 
fluctuating electronic structure along a molecular dynam- 
ics (MD) trajectory onto a low-dimensional model Hamil- 




E [eV] 

FIG. 1: Typical time-averaged transmission function {T{E)}^ 
for pG-pC and pA-pT wires containing seven base pairs each. 
The time series length was 100 ps and the average was per- 
formed every 5 ps. Inset; A snapshot of a pG-pC oligomer 
with 11 base pairs in a solvent drawn from the MD simulation. 
To avoid spurious boundary effects, only the innermost seven 
base pairs are used. Via the fragment orbital method, the 
electronic structure at each snapshot is mapped onto an effec- 
tive low dimensional model Hamiltonian with time-dependent 
electronic parameters. 
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FIG. 2: Coherence parameter Cn{E) which provides a qual- 
itative measure for the importance of dynamical effects; 
Cn{E) <^ 1 hints at the dominant role of fluctuations in de- 
termining the transport properties. Notice, however, that this 
parameter can only give reliable information when {T{E))^ 
0, i.e., within the energy window belonging to the spectral 
support of {T{E))^. The inset shows the energy averaged 
{Cn[E)) ^ as a function of the number of base pairs in the 
DNA chain. 



tonian. Charge transport will be studied for an effective 
model describing the coupling of the electronic system 
to a bosonic bath which comprises internal vibrations 
and solvent effects. The bath thus encodes dynamical 
information drawn from the MD simulations. Our treat- 
ment allows (i) the determination of electronic coupling 
parameters under realistic conditions, and (ii) the calcu- 
lation of the bath spectral density from time series gen- 
erated during the MD simulation. In this way, it does 
not contain any free parameters describing the molec- 
ular electronic structure or the coupling to the struc- 
tural fluctuations. We show, by considering as paradig- 
matic cases pG-pC and pA-pT oligomers, that the elec- 
trical transport properties in such bio-molecular systems 
are strongly dominated by the conformational dynam- 
ics. Despite some limitations of the model Hamiltonian 
approach related to the used approximations (see later), 
we nevertheless stress that the range of applicability of 
our method is not limited to DNA molecular wires; in- 
deed, it provides a solid basis for the parameter-free in- 
clusion of dynamical effects in a model-based treatment 
of quantum transport as well as for a multi-scale ap- 
proach to the description of the electronic properties 
of bio-molecules and their response to external fields. 
Our approach exploits a fragment orbital description [35 1 
of the bio-molecules, which allows for an efficient and 
well controlled coarse graining of the electronic structure 
problem [sil . Is^]- A combination of quantum mechan- 
ics/molecular mechanics methods (to describe solvent 
effects), MD simulations, and a parametrized density- 
functional tight-binding methodology [sl] has been used 
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FIG. 3; The time-averaged current {I{V))^ for a pGr-pCy 
oligomer. Results are shown for three different choices of the 
time frame tw =5, 20, and 50 ps (for a definition see the 
main text). Inset: Auto-correlation function of the onsite 
energy fluctuations (ei(f)ei(O)) and nearest-neighbor cross- 
correlation {ei(t)ei+i(0)). 



to extract the relevant electronic information in the form 
of time series. This leads to a time-dependent Hamilto- 
nian: H = V . e,{t)d]d, -f Y., V.-j+ii^) + l^'^')' 

where both ej(t) and Vj^j-\.i{t) are random variables as 
a function of the simulation time. These parameters 
describe, respectively, the effective ionization energy of 
a single base-pair —which defines a fragment in our 
calculations— and the coupling between nearest-neighbor 
fragments. We have approached the transport problem 
from two complementary perspectives. 

Time averaging and dynamical fluctuations — In Fig. 1 
we show the time-averaged transmission function {T{E))^ 
for p G-p C and pA-pT wires containing seven base 
pairs [39|. These calculations have been carried out for 
a Tmd = 30 ns long MD simulation with a time step 
of 1 fs. The first point to note is the apparently higher 
transmission of pA-pT compared to that of pG-pC. This 
is just the opposite of what a purely static calculation 
would yield. This fact represents a first hint at the im- 
portance of dynamical effects in determining charge prop- 
agation. The fragmented structure of the spectrum is 
simply mirroring the broad distribution of onsite ener- 
gies induced by the dynamical disorder. We have fur- 
ther defined a coherence parameter (CP) for a given 
chain length N as Cn{E) = (1 + ariE)/ {T{E)f^)-\ 
which can provide a quantitative measure for the role 
of structural fluctuations: Cn{E) <C 1 can indicate 
the dominance of the conformational dynamics. Hereby, 
aT{E) = {[T{E) - {T[E))^Y)^ is the mean-square devi- 
ation of the transmission and the brackets always indi- 
cate time averaging. Obviously, the CP has only a clear 
meaning in the spectral region where the transmission 
function itself is non-negligible (spectral support). The 



CP is shown in Fig 2 for seven base pairs of pA-pT and 
pG-pC, from where it is seen that (i) transport is dom- 
inated by the bio-molecular dynamics, Cn{E) <C 1, and 
(ii) the CP for pG-pC is roughly one order of magnitude 
smaller (within the spectral support domain) than for 
pA-pT, reflecting the fact that the latter system seems 
to be less affected by dynamical disorder. The inset 
of the same figure displays the energy-averaged CP for 
four different numbers of base pairs. Longer chains are 
clearly more affected by dynamical disorder than shorter 
chains, independent of the base sequence. In a second 
step, we have investigated the dependence of the time- 
average current {I(V))f. on the averaging procedure, i.e., 
calculating a set of partial currents Ii{V,tw) obtained 
upon averaging of the electronic parameters over time 
windows of length tw = ndSt along the time series, where 
(5i = 1 ps is the time step at which molecular confor- 
mations were extracted along the MD trajectory. The 
index £ ~ 1, • • ■ , int [Tmd/tw] = L labels the number 
of time frames once Ud has been fixed. The total cur- 
rent is thus given by {I{V))^ = {l/L)J2eIe{V,Tw). The 
different sizes of the time windows (different values of 
Tw) are mirroring in a phenomenological way differences 
in electronic time scales (an information not provided 
by the MD simulations); a charge will explore different 
fiuctuating environments in dependence of tw and thus 
the total current must be affected by this fact. Hence 
e.g., Tw <C with uj~^ being some typical time scale 
for dynamic fluctuations, would correspond to the non- 
adiabatic limit where a time-averaged atomic frame is 
felt by the charge, while the opposite limit Tyv ^ iiJ~^ 
defines the adiabatic regime, where instantaneous atomic 
configurations are "seen" . Of course, this provides only a 
qualitative picture, since the DNA structural fiuctuations 
involve many different time scales making the effective 
interaction of a charge with different degrees of freedom 
very complex. In Fig. 3, we show the time-averaged cur- 
rent for a fixed number of base pairs and three different 
values of the time frame: tw =5, 20, and 50 ps. We see 
a slow increase of the current with increasing tw, since 
a moving charge will effectively sample an increasingly 
larger number of realizations of ej{t) and Vjj+i{t). We 
remark that the current calculation using a Landaucr-likc 
expression is only meaningful! near the adiabatic limit; 
in the strongly non-adiabtic regime golden-rule like ex- 
pressions should be used [40| . 



Effective charge-bath model — As a complementary way 
to deal with charge transport —allowing for a flexible 
treatment of different transport mechanisms while still 
relying on a realistic description of the bio-molecular 
dynamics— an effective model has been formulated de- 
scribing the electronic system coupled to a fiuctuating 
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FIG. 4: I-V characteristics of pG-pC (upper panel) and 
pA-pT (lower panel) oligomers with three different num- 
bers of base-pairs A'^ = 7, 11, 15 calculated with the effec- 
tive charge-bath model Hamiltonian, see Eq. (1). In all cases 
the electrode-molecule coupling parameter was chosen P = 5 
meV. [4l|] Though for the shortest oligomer pA-pT has a 
larger current than pG-pC, this behaviour is weakened with 
increasing length. Possible reasons for this behavior are men- 
tioned in the text. 



environment (bosonic bath): 

+ J2x,jd]d,{B, + Bl) +J2^sBlB,. (1) 
Here, the time averages of the electronic parameters 



i^j)t ' (^J+i)t have been splitted off, e.g., ej{t) — (e 
6ej{t). Some approximations are involved by the formu- 
lation of this model: (i) only local energy fluctuations 
are considered and included in the bath (last term of 
Eq. (1)). (ii) Xsj (s numbers the bath modes) depends in 
general on the site j. This will be partially taken into ac- 
count and is reflected in a renormalization of the average 
hopping {Vjj+i)^; (iii) no fluctuations in the coupling 
parameters Vjj^i(t) are considered. The bath will be 
characterized by a site-averaged spectral density J{uj) = 
{6e^{0)) (2/7rft) tanh (fiw/fceT) dt cos (ujt) C'{t), with 
C{t) = il/N)J2j {Sej(t)Sej{0)) being the (site-averaged) 
autocorrelation function of the onsite energy fluctuations. 
Using the model of Eq. ([1]), the electrical current through 
pG-pC and pA-pT oligomers containing 7V=7, 11, and 15 
base pairs was computed [4l| . The Fermi level was fixed 
in each case outside the spectral support of T{E) to ob- 
tain a zero-current gap, which is thus arbitrary in these 
calculations. In Fig. 4, where the I-V characteristics of 
the different sequences and lengths are shown, we observe 
that apart from the shortest {N = 7) oligomer, the cur- 
rent for pG-pC is somewhat larger than for pA-pT. These 
features are possibly associated to two factors. Firstly, 
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the neglect of non-local onsite energy fluctuations: fluc- 
tuations between neighboring sites decay faster on sub-ps 
time scales, see the inset of Fig. 4, but can nevertheless 
induce non- vanishing correlations over few base pairs (37j 
thus modifying the electrical response of the system. Sec- 
ondly, the use of an averaged spectral density effectively 
makes the coupling of all electronic sites to the bath very 
similar. Since fluctuations become more important with 
increasing length, this approximation may become prob- 
lematic. In spite of these limitations, Eq. ^ provides a 
reasonable starting point to bridge MD simulations with 
charge transport models and, more important, offers the 
possibility of systematically improving the model Hamil- 
tonian approach. 

In conclusion, we have presented a combined molecular 
dynamics / model Hamiltonian approach which allows for 
a very flexible treatment of charge transport through bio- 
molecular systems taking into account dynamical disor- 
der. The method can allow, as illustrated in the special 
case of DNA wires, the straightforward study of the base- 
sequence and length dependence of the electrical response 
of such systems. The results presented here strongly sup- 
port the view that charge transport through DNA wires is 
dominated by conformational fluctuations. In this sense, 
transport approaches based on band-like coherent trans- 
port or on purely static structures can not yield a real- 
istic description of charge motion in such highly dynam- 
ical systems. Finally, we would like to stress that our 
approach can be applied as well to investigate the in- 
terplay of charge transport / transfer and conformational 
dynamics in other complex bio-molecular systems. This 
flexibility relies on the fact that the degree of coarse- 
graining —leading to the formulation of effective tight- 
binding models— can be "tuned" by an appropriate re- 
definition of the fragment orbitals while still retaining the 
relevant dynamical information. 
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